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ABSTRACT 

We report upper limits on the masses of black holes that can be present in the centers of sixteen nearby galaxy 
bulges. These limits Mg' H for our statistically complete sample were derived from the modeling of the central 
emission-line widths ([N II] or [S II]), observed over a 0"25 x 0"2 (R < 9 pc) aperture. The experiment has a 
mean detection sensitivity of ~ 3.9 x 1O 6 M . For three sample members with direct determinations of Mbh our 
upper limits agree within the uncertainties, while in general our upper limits are found to be close to the masses 
measured in other bulges with global properties similar to ours. Remarkably, our limits lie quite closely to the 
recently derived Mbh _ relation. These results support a picture wherein the black-hole mass and overall galaxy 
structure are closely linked, as galaxies with exceptionally high Mbh at a given o± apparently are rare. 

Subject headings: galaxies: bulges — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: spiral 



1. INTRODUCTION 

The last few years have seen great progress in studying the 
dark mass concentrations at the centers of "ordinary" quies- 
cent galaxies, showing that they are very common and demon- 
strating in some cases that they must be supermassive black 
holes (SMBHs) by ruling out all astrophysically viable alter- 
natives. Indeed, a picture is emerging in which SMBHs in the 
range of 10 6 to 10 9 M Q are an integral part of galaxy formation 
(e.g. , Kauffmann & Haehnelt 2000). Two examples in the lo- 
cal universe stand out with particularly convincing evidence as 
SMBHs: the Milky Way and NGC 4258. At the Galactic Cen- 
ter, direct observations of individual stars (Genzel et al. 1997; 
Eckart & Genzel 1997; Ghez et al. 1998) and a stream of ion- 
ized gas (Herbst et al. 1993) orbiting Sgr A* show that all the 
dynamically relevant mass inside ~ 1 pc, 2.6 x 1O 6 M , is con- 
centrated with a density of p > 10 12 M Q pc" 3 . In NGC 4258, 
a disk of masing molecular gas is orbiting the center with 
a Keplerian rotation curve as traced by H2O maser emission 
(Miyoshi et al. 1995); models imply thatM BH ~ 3.6 x 1O 7 M 
and p > 4 x 10 9 M Q pc" 3 . Alternatives to SMBHs, such as clus- 
ters of brown dwarfs or stellar remnants, can be ruled out in 
these two cases (e.g. Maoz 1995, 1998). 

A number of techniques have matured that have demon- 
strated the presence of a central dark mass concentration in an 
ever growing number of nearby galaxy nuclei, with mass es- 
timates accurate to a factor of ^2 and concentration limits of 
p > 10 6 ~ 8 Mq pc" 3 . Simple analogy with the exemplary cases 
where the SMBH presence is all but proven, and the connec- 
tion to active galactic nuclei (AGN) activity at various inten- 
sity levels, suggest strongly that these central dark masses are 
SMBHs as well. The most widely used technique is stellar dy- 



namical modeling (e.g., Dressier & Richstone 1988; Kormendy 
et al. 1996, 1997; van der Marel et al. 1997; Cretton & van den 
Bosch 1999; Gebhardt et al. 2000a), which has provided mass 
estimates for over two dozen nuclei, mostly in massive, early- 
type galaxies. Modeling the kinematics of ionized gas has pro- 
duced a number of additional mass measurements (e.g., Harms 
et al. 1994; Ferrarese, Ford, & Jaffe 1996; Bower et al. 1998; 
Verdoes Kleijn et al. 2000; Sarzi et al. 2001, hereafter S01; 
Barth et al. 2001a). Finally, application of results from rever- 
beration mapping of active galactic nuclei have yielded central 
virial mass estimates for Seyfert galaxies (Ho 1999; Wandel, 
Peterson, & Malkan 1999) and QSOs (Kaspi et al. 2000) that 
seem to be robust (Gebhardt et al. 2000b; Ferrarese et al. 2001). 

Taken together, these results show that Mbh is correlated both 
with the stellar luminosity Lbuige (Kormendy & Richstone 1995; 
Magorrian et al. 1998; Ho 1999; Kormendy et al. 2001) and, 
more tightly, with the velocity dispersion of the bulge o± (Geb- 
hardt et al. 2000a; Ferrarese & Merritt 2000). By contrast, M BH 
is unrelated to the properties of galaxy disks (Kormendy et al. 
2001). The growth of SMBHs appears to be closely linked with 
the formation of bulges. However, the actual slope and scat- 
ter of the Mbh _ ^buige and Mbh _ &± relations are still under 
debate. It is also important to remember that our knowledge 
about SMBHs is very uneven across the Hubble sequence of 
galaxies. The existing samples are preferentially weighted to- 
ward early-type galaxies with very massive black holes. From 
an observational point of view, there is a pressing need to ac- 
quire better Mbh statistics for spiral galaxies. 

Motivated by the recent progress and the emerging correla- 
tions, but also by the desire to improve the black hole census 
in spirals, we derive mass constraints on SMBHs potentially 
present in the bulges of sixteen nearby disk galaxies. As we 
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will show, these constraints constitute significant progress, both 
in terms of the number of target galaxies and in terms of broad- 
ening the range of parent galaxies with significant constraints. 

We draw on spectra of nearby nuclei obtained with the Space 
Telescope Imaging Spectrograph (STIS) onboard the Hubble 
Space Telescope (HST), taken as part of the Spectroscopic Sur- 
vey of Nearby Galaxy Nuclei (SUNNS) project (Shields et al. 
2000; Ho et al. 2000; S01; Rix et al. 2001). Only four of 
our original twenty-four target galaxies showed extended line 
emission with symmetric kinematics and hence were suited for 
a direct M B h determination. These cases have been modeled by 
S01 and yielded M Bli estimates. Here we analyze data for six- 
teen galaxies from SUNNS with central [N II] AA6548, 6583 
and [S II] AA6716, 6731 line emission from ionized gas. In 
most cases the central gas emission is spatially unresolved, or 
only marginally so, and has line widths <~ 10 2 kms" 1 . 

Spatially unresolved lines do not permit precise mass esti- 
mates, but potentially do allow us to derive useful upper limits 
on Mbh- This is because the line emission in our central aper- 
ture must arise from gas at a distance from the center that is, at 
most, equal to the physical dimension of the region subtended 
by the central aperture itself, for our sample typically ~ 9 pc. 
If the gas motions are orbital, all velocities and hence the in- 
tegrated line width will scale as \/Mbh- Note that the result- 
ing limit on Mbh scales linearly with the central aperture size, 
affording HST an order of magnitude gain over ground-based 
observations (e.g., Salucci et al. 2000), and making the derived 
limits astrophysically interesting. 

At parsec-scale distances from galactic centers the gas is sub- 
ject not only to gravitational forces but also to gas pressure 
and magnetic forces (e.g., for the Milky Way, Timmermann et 
al. 1996; Yusef-Zadeh, Roberts, & Wardle 1997). In general, 
these other effects cause additional line broadening. 

The paper is organized as follows. In §2 we present the spec- 
troscopic and photometric STIS observations, and in §3 we de- 
scribe our modeling of the ionized gas kinematics. In §4 and 
§5, we present our results and draw our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Galaxy Sample 

All sample galaxies, mostly early-type disk galaxies (SO - 
Sb), were observed with STIS as part of the SUNNS project; 
the full details of this program will be reported elsewhere (Rix 
et al. 2001). The SUNNS galaxies are drawn from the Palomar 
spectroscopic survey of nearby galaxies (Filippenko & Sargent 
1985; Ho, Filippenko, & Sargent 1995, 1997) and include all 
SO to Sb galaxies within 17 Mpc known to have line emis- 
sion (> 10 15 ergs s" 1 cm" 2 ) within a 2" x 4" aperture. The 
present sample constitutes the subset of SUNNS objects with 
sufficient central line flux in Ha or [N II] to provide adequate 
signal-to-noise (S/N > 10). With the experimental setup de- 
scribed below, the sample is effectively defined by galaxies that 
have Ha or [N II] line fluxes > 10~ 14 ergs s" 1 cm" 2 within a 
0."25 x 0"2 central aperture, which correspond to line lumi- 
nosities of > 3 x 10 38 ergs s" 1 at the mean sample distance of 
15 Mpc. Our sample is statistically well defined and complete 
in the sense that the original SUNNS sample was a volume- 
limited sample selected by emission-line flux within a (much 
larger) 2" x 4" aperture. The basic parameters of the target 
galaxies are given in Table 1 . 

1 1 IRAF is distributed by the National Optical Astronomical Observatories, which 



2.2. Observations 

HST observations were acquired for all objects in SUNNS 
during 1998 and 1999. We placed the 0"2 x 52" slit across 
each nucleus along an operationally determined position angle, 
which is effectively random with respect to the galaxy orien- 
tation. After initial 20-s acquisition exposures with the opti- 
cal long-pass filter (roughly equivalent to R), from which we 
derive surface photometry of the central regions, three expo- 
sures totaling approximately 45 minutes were obtained with the 
G750M grating; this resulted in spectra that cover 6300 A to 
6850 A with a full-width at half maximum resolution for ex- 
tended sources of 1 .6 A. 

For 9 of the 16 galaxies in the present sample, the tele- 
scope was offset by 0."05 (~1 pixel) along the slit between 
repeated exposures to aid in the removal of hot pixels and cos- 
mic rays. The two-dimensional (2-D) spectra were bias- and 
dark-subtracted, flat-fielded, aligned, and combined into single 
frames. Cosmic rays and hot pixels remaining in the combined 
2-D spectra were cleaned following the recipe of Rudnick, Rix, 
& Kennicutt (2000). The 2-D spectra were then corrected for 
geometrical distortion and then wavelength and flux calibrated 
with standard stsdas procedures within iraf 11 

To represent the generic "nuclear spectrum" of each galaxy, 
we extracted aperture spectra five pixels wide (<~ 0"25), cen- 
tered on the brightest part of the continuum. The acquisition 
images indicate that the uncertainties in the galaxy center due 
to dust are <; 0.25 pix <~ 0."012 (see also S01). In essence, 
therefore, the extracted spectra represent the average central 
emission, convolved with the STIS spatial point-spread func- 
tion (PSF) and sampled over an aperture of 0"25 x 0"2, or 1 8 pc 
x 14 pc for the mean sample distance of 15 Mpc. 

For three of our galaxies, central stellar velocity disper- 
sions were either not available in the literature (NGC 3992 
and NGC 4800) or quite uncertain (NGC 3982; Nelson & 
Whittle 1995). Therefore, we obtained new spectroscopic data 
for these objects on two observing runs: with the Boiler & 
Chivens spectrograph at the Bok 90-inch telescope in May 2000 
for NGC 3992, and with the Double Spectrograph (Oke & 
Gunn 1982) at the Palomar 200-inch telescope in June 2001 
for NGC 3982 and NGC 4800. At Kitt Peak we used the 600 
grooves mm" 1 grating to cover 3600-5700 A with a pixel scale 
on the CCD of 1.86 A pix" 1 , while for the Palomar run we 
observed the Ca infrared triplet using the 1200 grooves mm" 1 
grating on the red side of the Double Spectrograph, with a pixel 
scale of 0.63 A pix" 1 . Spectra were extracted for apertures of 
3 "3 x 2."5 and 3 "7 x 2."0 for the Bok and Palomar spectra, 
respectively. Total exposure times were 40, 30, and 60 min- 
utes for NGC 3982, NGC 3992, and NGC 4800, respectively. 
The stellar velocity dispersions were measured following the 
method of Rix et al. (1995), and the obtained values are re- 
ported in Table 1 . 

2.3. Central Emission-Line Widths and Flux Profiles 

To quantify the emission-line velocity widths in these nuclear 
spectra, we simultaneously fit Gaussians of single width a cen 
to the [N II] AA6548, 6583 and [S II] AA6716, 6731 emission- 
line doublets, using the iraf task specfit. The profiles of the 
[N II] and [S II] lines were found to be roughly Gaussian. We 
restricted ourselves to line widths from forbidden transitions 
to side-step the impact of a possible broad (broad-line region) 

operated by AURA, Inc. under contract to the NSF. 
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component arising from radii much smaller than our observa- 
tional aperture. However, in all objects where only a narrow 
Ha emission-line component was present, we also included Ha 
in the fit. In the few cases with prominent, very broad Ha 
lines (e.g. , NGC 4203, Shields et al. 2000; NGC 4450, Ho 
et al. 2000), particular care was taken to minimize the impact 
of the very broad lines on the estimate of a cen in the adjacent 
[N II] lines. In virtually all objects the S/N is in excess of 10, 
and hence the formal errors in the estimated line width are neg- 
ligible for the subsequent analysis. The instrumental line width 
derived from comparison lamps is a inst s=a 32 kms" 1 and was 
subtracted from the raw measurement of a cen in quadrature; for 
all but two objects these line width corrections were negligible, 
implying that the intrinsic widths were well resolved. This cor- 
rection also spares us from accounting in the following model- 
ing of the line width for the broadening due to the instrumental 
line spread function. The resulting values for a cen are listed in 
Table 1 . Their characteristic errors, including the correction for 
the instrumental line width, are less than 10 kms" 1 . 

As we will detail below, any information on the gas spatial 
flux distribution on scales < 0."25 provides a valuable con- 
straint for the central line width modeling procedure. There- 
fore, we also obtained radial profiles along the slit direction for 
the ionized gas flux, by fitting a Gaussian to the [N II] A6583 
line-flux profile along the spatial direction on the 2-D spectra. 
We chose [N II] because among our sample galaxies this is al- 
most always the brightest line, and because it is less likely to be 
affected by underlying absorption features in the stellar contin- 
uum than Ha. The [N II] emission-line flux profiles are shown 
in Figure 1 . 



3. MODELING THE CENTRAL LINE WIDTH 

3.1. Basic Concept 

We are now faced with converting the observed central line 
widths into estimates for the central black-hole mass. 

To start, we assume that the ionized gas motion is dominated 
solely by gravity. In this case the central line width depends on: 
(a) the total gravitational potential of the putative SMBH and 
of the surrounding stars; (b) the spatial emissivity distribution 
(e.g. , that of a disk inclined at cos/); and (c) the "kinematic 
behaviour" of the ionized gas, for example "dynamically cold" 
gas moving on circular orbits or hotter gas with hydrostatic sup- 
port. The lack of spatially resolved information on the gas flux 
distribution within the central 0."25 x 0"2 aperture means that 
we can only derive upper limits to Mbh; if the emission-line 
flux within the aperture arose from R <C Rapemm « 0."1, ar- 
bitrarily small values of M B h could explain the observed line 
width. If the gas motion is also affected by non-gravitational 
forces, such as outflows, magnetic fields, or supernova winds, 
this would broaden the integrated line velocity width addition- 
ally, and hence lower the required black-hole mass needed to 
explain a given a cen . By ignoring non-gravitational forces we 
are therefore conservative in estimating upper limits for the cen- 
tral black-hole mass. The absence of constraints on the impor- 
tance of non-gravitational forces constitutes a second reason 
(besides the lack of spatially resolved information on the gas 
flux) why it is not possible to secure the presence of a SMBH in 
our sample galaxies, as hypothetically the observed line widths 
could be entirely explained by non-gravitational effects. 

If the functional form of the potential well is fixed, then the 
central line width will scale with the potential for any given 



choice of the emissivity distribution and for the gas kinematical 
behaviour. In the simplest case of a purely Keplerian potential 
induced by a SMBH, $ B h, the expected central line width will 
scale as the square root of the black-hole mass. As the circular 
velocity at any given reference radius R re f, v c (R re f), will scale 
in the same way, the ratio between o cm and v c (R re f) is inde- 
pendent of black-hole mass. The task at hand is therefore to 
derive a plausible range of values for this ratio by varying the 
spatial emissivity distribution, and then to obtain a mass range 
for the putative SMBH from the observed central line width via 
v 2 c (R ref ) -» M BH = v 2 c (R ref )R ref /G. 

The same would hold for a sequence of purely stellar poten- 
tials derived from the luminosity density with differing mass- 
to-light ratios T. When both the stellar and the SMBH contri- 
bution to the gravitational potential are considered, the shape of 
the rotation curve, and hence a cen /v c (R re f), will depend on the 
relative weight of M Bll and T. 

In all cases we will proceed through the following steps in or- 
der to make a prediction for the gas velocity dispersion within 
the central aperture: 

• Specify the spatial gas emissivity distribution and the 
gravitational potential and choose the kinematic be- 
haviour of the gas. 

• Construct 2-D maps for the moments of the line-of-sight 
velocity distribution (LOSVD) at any position (x,y) on 
the sky 

Sv*(jc,y) = J LOSVD (x,y,v,)v k z dv z (k = 0,1,2), (1) 

as they would appear without the limitations of the spa- 
tial resolution; the first moment, for instance, is the mean 
gas velocity. 

• Convolve each of the 2-D £v* maps with the STIS PSF. 

• Sample the convolved Y,v k conv 2-D maps over the desired 
aperture to obtain the PSF-convolved, aperture-averaged 
LOSVD velocity moments Y>v k comKap , which are directly 
comparable to the observables. 

• In particular, compute the ionized gas flux f ap , 
the projected mean streaming velocity v ap , and the 
velocity dispersion a ap within the desired aperture 

th rough f ap = £y° co „ w , v ap = Sv 1 conv .ap / 'f ap , and <J ap = 

\J Ev 2 comvip /f ap -Va P , respectively. The last quantity, 
c a p, can be compared with the measured velocity width 

&cen- 

In what follows, we will first derive upper limits on Mbh as- 
suming that the gas is moving on circular orbits in a coplanar, 
randomly oriented disk within a Keplerian potential (§3.2); then 
we consider the impact of the stellar potential on this disk mod- 
eling (§3.3); finally, we will explore a seemingly very different 
situation for the kinematical behaviour of the gas, that of hy- 
drostatic equilibrium (§3.4), to demonstrate that our results are 
robust with respect to the underlying model assumptions. We 
anticipate that the most conservative upper limits on M B h are 
derived from the first approach, where the impact of the stellar 
potential is neglected. 
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FIG. 1. — Observed emission-line flux profiles along the direction of the STIS slit. In each panel the filled symbols represent the observed [N II] fluxes at each 
CCD row along the slit (corresponding to a synthesized aperture of 0"05 X 0"2), while the dashed line shows the fit to the central five points as described in §3.2. 
For comparison the solid line represents the stellar continuum profile, extracted from the region free of emission lines between 6536 and 6542 A. Both emission-line 
and continuum flux profiles have been rescaled so that their peak values equal unity. 



3.2. The Keplerian Disk Modeling 

We start with the simple and plausible assumption that the 
ionized gas moves on circular orbits at the local circular ve- 
locity, which in turn is dictated solely by the gravitational in- 
fluence of the putative SMBH, v 2 c {R) = GM BH /R. We further 
assume that the gas resides in a coplanar disk of unknown in- 
clination with an intrinsically axisymmetric emissivity distribu- 
tion, £(/?), centered on the stellar nucleus. Our best guess for 
£(/?) is derived from the data themselves. In this Keplerian disk 
the LOSVD at each position (x,y) on the sky plane is just 

LOSVD(x,y, v z ) = ^ P roj(x,y)S[v CtProj (x,y)-v z ] (2) 

and its T,v k velocity moments are simply given by Y, pro j(x,y), 
Z pm j(x,y)v CtPn)j (x,y), and 5Xx,y)v% iPmj (x,y), respectively, 
where T, pro j and v CtPro j are the projected gas surface brightness 
and circular velocity. To deal with the central circular velocity 
singularity, we neglected the contribution of the central point 
of the 2-D maps for the Sv 4 velocity moments, and we refined 
their grid sizes until no further substantial increase in the pre- 
dicted line widths was found. The adopted grid size in our mod- 
els corresponds to 0"005, or one tenth of a pixel. 

The geometry of the projected velocity field will depend on 
the disk orientation, specified by its inclination i with respect 
to the sky plane and its major axis position angle <f> with re- 
spect to the slit direction. We have no information on the gas 
disk orientation within our central aperture, since the dust-lane 
morphology we employed in S01 cannot be used to provide a 



constraint on such small scales. Therefore, we need to explore 
all possible disk orientations to derive the probability distribu- 
tion for cr ap /v c (R re f), the ratio between the predicted central 
velocity dispersion and the circular velocity at the reference ra- 
dius. We adopt R re f = 0."125, corresponding to the distance of 
the central aperture edge from the center along the slit direc- 
tion. We cover the possible disk orientations by constructing a 
grid of models with equally spaced cos i and <f). 

For the intrinsic radial surface brightness profile of the gas 
we assumed a Gaussian 

nR) = a f i ux e- Rl l 2(j2 f"' % (3) 

and we derived a// M and afi m by matching the observed emis- 
sion flux profiles along the slit (see §2.3). This match again in- 
volves convolving the intrinsic with the STIS PSF, which 
we parameterized as a sum of Gaussian components (see S01). 
The choice of a Gaussian for the intrinsic surface brightness 
distribution was a matter of convenience, for the convolution 
process in this case is simply analytical. Intrinsically more con- 
centrated profiles, such as exponential ones, would also repro- 
duce the data once convolved with the STIS PSF, and would 
lead to tighter upper limits on M B h, which make our choice 
more conservative. We fit only the central five flux pixels for 
each galaxy, corresponding to the same region subtended by 
our central aperture. 

Figure 1 displays our best fits of £(/?), and illustrates that our 
model (Eq. 3) matches the data well within < 0."125 in all 
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FIG. 2. — Predicted ratio of the gas central velocity dispersion to the circu- 
lar velocity at the reference radius, as obtained by the Keplerian disk modeling 
(§3.2) for different disk inclinations, position angles <f>, and intrinsic flux dis- 
tributions. 

cases. 12 Table 2 lists the best-fitting Ofi ux values. The actual 
position of the data points with respect to the center of the fit- 
ting function can be explained entirely by a small displacement 
(<C Rref) °f the slit center from that of the gaseous disk, without 
violating our assumption of an axisymmetric emissivity distri- 
bution. Including this small offset in the modeling produces 
only negligible variations in the predicted central u ap . In this 
case the predicted mean velocities v ap , although non-zero, are 
very small and consequently the velocity dispersions a ap are 
almost equal to the ones obtained with perfectly centered aper- 
tures (when v ap = 0). 

For comparison, we also show the stellar surface brightness 
profiles in Figure 1, derived from the stellar continuum in the 
spectra; this comparison justifies our assumption that the gas 
and the stellar distribution are concentric. Indeed, an indepen- 
dent fit to the stellar profiles with the same functional profile 
adopted for the gas ones led to a mean offset along the slit di- 
rection of only 0.08 ± 0.17 pixels, consistent with no offset at 
all. 

As an intermediate result of our modeling, we show in Fig- 
ure 2 the <J av /v c {R re f) ratios obtained in a Keplerian potential 
for different disk orientations and a range of typical values for 
(jfi U x\ the predicted <j ap does not depend on the total flux, or 
a fiux. The predicted central line width obviously increases from 
face-on to edge-on systems. At a given disk orientation, mod- 
els with intrinsically more concentrated gas emissivity always 
have a larger line width than those with more extended flux dis- 
tributions because the gas resides at smaller radii. Since the 
central 0."25 x 0"2 aperture is nearly square, the impact of the 
position angle parameter <fr on the final confidence limits for the 
0a/>/v e (^ref) rat i° i s negligibly small. 

The flux distributions in our sample galaxies are concentrated 
enough that the predicted line widths are monotonically 
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decreasing with increasing cos i. Since randomly oriented disks 
have uniformly distributed cos i, we can use Figure 2 to derive 
the median values and 68% upper and lower confidence lim- 
its for the Mbh, respectively, by simply taking the values of 
&a P /v c (Rref) for the models with cos; = 0.5, 0.84, and 0.16. 

There is a minor practical complication in this modeling: 
while changing the orientation of a disk with a fixed intrinsic 
flux distribution, the predicted flux profile along each of the five 
0."05 x 0."2 apertures within the central 0."25 x 0."2 aperture 
changes, eventually becoming inconsistent with the observed 
one. Hence, at any given disk orientation we must readjust the 
intrinsic flux concentration <jfi ux in order to match the central 
five flux data points. Such a correction is particularly important 
for highly inclined disks, which have very different flux profiles 
when considered at different position angles. Since our M B h 
upper limits are derived for nearly face-on orientations, these 
corrections will not strongly affect our results. Furthermore, 
tests have shown that the induced scatter in the a ap jv c {R re f) 
ratio at any given inclination is considerably smaller than the 
face-on to edge-on variation. 

For simplicity, we only used a statistical correction for this 
effect. For a set of intrinsic values a/iuxjn, we collected all the 
central flux profiles predicted for a uniform grid in cos; and 
<f>. Then we treated each of these profiles as observed ones and 
matched each of them with a PSF-convolved Gaussian profile 
to get a distribution of afi ux . oul values and a median (<Jfi ux . out ) 
value for each a f i ux . in . By comparing the median (a f i ux . out ) 
values with the corresponding Ofi m j„ values (Fig. 3), we can 
correct for each galaxy our initial guess of the intrinsic a/i ux , 
derived by matching the observed central flux profiles of 



0.1 

(arcsec) 



0.15 



12 All [N II] flux profiles of Figure 1 were symmetrical even outside the central aperture region, with the noticeable exception of NGC 4501 and NGC 4698. This 
last Sa galaxy shows the presence of a stellar (Bertola et al. 1999) and gaseous (Bertola & Corsini 1999) core with an angular momentum perpendicular to that of the 
main galactic disk. This core can be identified as a disk from HST imaging (Scarlata et al. 2001). Consistent with these findings, a recent accretion event may explain 
why NGC 4698 is one of the two galaxies, among our sample of 16, to exhibit a strongly asymmetric gas distribution within 0"3 (±6 rows) where, assuming that the 
nuclear regions of our galaxies are dominated by SMBHs with masses consistent with the Mbh -o* relation, the dynamical timescale ranges from 0.5 to 12.0 Myr. 
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FIG. 4. — Physical dimension of the region subtended by our central 
()"25 X 0'.'2 aperture (filled circles) and flux distribution extent of our sample 
galaxies (open stars) compared to the "sphere of influence" radius of SMBHs 
consistent with the Mbh _ relation (solid line). 

Figure 1. For a given galaxy the corrected flux concentra- 
tion to be input into the models is characterized by the (Jfiuxjn 
value that, according to the previous scheme, lead to a median 
(<Jfiux,out) equal to the a flux of that galaxy. This cr /te , ( „ value 
describes the intrinsic flux distribution that, when considering 
different disk orientations, leads to the predicted flux profiles 
that are the most consistent with the observed one for the given 
galaxy. 

We call these particular Ufi^n for each galaxy the corrected 
Pfiux.corr values, which we list in Table 2 along with the (+lcr) 
upper limits on M B h obtained adopting them. 

As a final remark we notice that at the mean distance of our 
sample galaxies (^15 Mpc) and considering the typical mass 
of their central SMBH's (~ 2.4 x 10 7 M Q , as predicted by the 
Mbh - 0* relation), the adopted intrinsic flux distributions are 
generally concentrated enough that a double-horned LOSVD 
should be expected, expecially for highly inclined disks. The 
fact that such a feature is not found in the observed emission- 
line profiles could represent evidence of an intrinsic turbulence 
in the gas. Indeed, when the full LOSVD is properly con- 
structed (by collecting at each velocity bin the total flux within 
the 0"25 x 0."2 aperture that arises from the corresponding iso- 
velocity slice of the flux distribution, once convolved with the 
STIS PSF), the double-horned shape disappears when an intrin- 
sic gas velocity dispersion is introduced. In the favorable case, 
within this context, of a nearly face-on disk with cos; = 0.84, a 
velocity dispersion of about <~ 30 kms" 1 would be required on 
average to smooth the double-horned shape, thus increasing the 
predicted line widths by 9% and decreasing the derived upper- 
limits on Mbh by 19%. Too many assumptions will be required 
to make this correction on the case by case basis, and would 
probably require direct fitting of the observed emission lines, 
as done by Barth et al. (2001b). By adopting our current ap- 
proach, our limits on M BH remain conservative upper-bounds. 

3.3. The Stellar Contribution 

We now proceed to evaluate the impact of the stellar po- 
tential on our modeling. We start by estimating for each 



galaxy the expected radius of the "sphere of influence" of the 
SMBH, r in fi = GM Bii /al, within which M B h dominates the dy- 
namics of a galaxy with a stellar velocity dispersion cr*. For 
this estimate we adopt the M BH - a+ relation as parameterized 
by Gebhardt et al. (2000), logM BH = 3.751ogo-*-0.55. Then, 
log r in f i = 1.75 log ov-2. 92, where M Bli is in units of M Q , a* in 
km s" 1 , and r in fi in pc. In Figure 4 we compare r in fi with the 
physical scale corresponding to the mean radius of our central 
aperture R aperture = ^(0"25 x 0"2)/tt w 0." 1 3, and find that for 
many of our galaxies r,„// < R aper ture, indicating that the stel- 
lar mass within R ap erture is comparable to or exceeds M Bli . 
Therefore, our analysis needs to account for the stellar mass, 
and we need to derive the stellar mass density profiles ^*(r), in 
particular for galaxies with smaller ov 

The relative importance of M* and M Bli on the central line 
width depends also on the spatial extent of the gas emissivity. In 
particular even when r,„/; < R aper ture, the presence of the SMBH 
can still be noticed from the observed a cen if Ofiux ~ Hnfi — 
that is, if most of the collected flux within R ape rture was emitted 
by gas moving on nearly Keplerian orbits. Further, the impact 
of on the predicted a ap also decreases with increasing flux 
concentration since in general the circular velocity curve due 
to the stars increases monotonically with the distance from the 
galactic center. Indeed the central slope of stellar densities pro- 
files can be represented by a power-law v+(r) ~ r~ a with a < 2 
(Gebhardt etal. 1996). 

Figure 4 shows ct/; m is in general smaller than R ape rture, an d 
we expect that the inclusion of the stellar mass in our modeling 
will cause only a modest correction in the black-hole masses 
inferred in §3.2. 

To quantify the stellar mass contribution in each galaxy, we 
derived the mass density profile i/*(r) by deprojecting the stellar 
surface brightness distribution £*(/?) obtained from the STIS 
acquisition image, assuming spherical symmetry and a constant 
mass-to-light ratio T. We applied the same multi-Gaussian al- 
gorithm adopted in S01 to circularly averaged Y,*(\fab) sur- 
face brightness profiles, extracted using the iraf task ellipse and 
color corrected into Johnson /?-band magnitudes using the iraf 
package synphot and assuming E-S0 galaxy templates. Gaus- 
sian components with a < 0.5 pixel were considered as unre- 
solved point sources and hence were excluded from the stel- 
lar mass budget. For simplicity we adopted for all galaxies 
T = 5 M Q /L Q , rescaled from van der Marel (1991) for H = 75 
km s" 1 Mpc" 1 , instead of deriving individual values for T by 
matching ground-based cr* measurements (see S01). 

Including the stellar potential in the modeling results in a 
27% reduction of the median black-hole masses (for cos ; = 0.5) 
needed to explain the observed central line widths. Nonethe- 
less, this effect is important for galaxies with small cr* values; 
the median M Bll decreased by 37% for the sample galaxies with 
the lowest a* (NGC 3982, NGC 4321, and NGC 4548), but only 
by 3% for the ones with the highest cr* (NGC 4143). 

The impact of the stellar mass on the upper limits of the 
black-hole mass, which are central to our analysis, is yet 
smaller, on average less then 12% (see Tab. 2). Indeed, similar 
to the purely Keplerian case, the predicted central line widths 
are always found to be increasing monotonically with disk in- 
clination, such that the derived +la upper limits on Mbh are 
obtained here also from models with nearly face-on disks. In 
this situation the circular velocities needed by the model to ex- 
plain the observed line widths always by far exceed the ones 
provided only by the stellar potential. 
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FIG. 5. — Mbh upper limits in the Mbh — c+ relation. The +lcr upper limits 
on Mbh (thick downward arrows ) obtained by the Keplerian disk modeling 
(§3.2) are compared with the range of black-hole masses expected as a func- 
tion of <r* from the Mbh — c* relation and its absolute scatter in Mbh (solid 
and dashed lines , from Gebhardt et al. 2000a). The open squares show the 
Mbh measurements obtained from extended kinematics in four of our sample 
galaxies by S01, while the open circles represent the sensitivity limit of our 
experiment for each sample galaxy. The thick horizontal error bars indicate 
how the mean error on it* for our sample would appear in the right and left 
hand of the plot. 

3.4. Gas in Hydrostatic Equilibrium 

So far we have assumed that the observed central line widths 
are due to pure orbital motion of gas in a disk around the galaxy 
center. But even if we ignore non-gravitational effects, we still 
need to investigate how much the derived Mbh upper limits de- 
pend on different choices for the kinematic behaviour of the 
ionized gas. In particular, it is conceivable that the gas is at least 
in part supported by gas pressure. In order to explore the im- 
pact of such pressure on our results we considered the extreme 
case of pure hydrostatic support. For any emissivity density 
profile p, we need only the second losvd velocity moment Sv 2 
without streaming motions. We obtain this moment by solv- 
ing the hydrostatic equilibrium equation d(pa 2 ) / pdr = —d§ jdr 
for the gas velocity dispersion cr(r), and then by integrating the 
luminosity-weighted a along the line of sight. 

Under these assumptions, we proceeded to match the ob- 
served <j cm for all the galaxies, assuming Gaussian profiles for 
their emissivity densities p, whose width has to match the ob- 
served central flux profiles. We thus obtained for each galaxy a 
value for the mass of the putative SMBH, in a purely Keplerian 
potential. 

From this exercise we found that the predicted central veloc- 
ity dispersions a ap correspond to values from the rotating disk 
model with disk inclinations always between 0.65 and 0.71. 
Consequently, the M B h values inferred by considering hydro- 
static equilibrium lie within the ±la confidence limits obtained 
from the rotating disk models, and all previous results hold. 

4. RESULTS 

We have explored the dynamical implications of the observed 
emission-line widths arising from the central <~ 10 pc of our 
sample galaxies, and we have demonstrated that the most con- 



servative lcr upper limits of Mbh are obtained assuming that 
the gas resides in a nearly face-on disk (z ~ 33°; cos; = 0.84) 
moving in circular orbits around a central SMBH. Fortunately, 
extremely face-on orientations are statistically rare. 

In Figure 5 we place our Mbh upper limits derived for the 
Keplerian case (§3.2) in the Mbh-ct* plane . For a comparison 
with the M B h - cr* relation (as parameterized by Gebhardt et 
al. 2000a) we need to scale each cr* to a e , the value that would 
have been observed within a circular aperture of R = R e . The 
proper computation of such a quantity, which should include 
the contribution from the rotation, was not possible for all the 
galaxies in our sample, since the necessary combination of sur- 
face brightness, velocity dispersion, and rotation velocity radial 
profiles was not available from the literature for all objects. The 
only choice in this case is to use the algorithm of J0rgensen, 
Franx, & Kjaergaard (1995), which was derived on the basis 
of kinematical data for 5 1 elliptical and lenticular galaxies and 
which also accounts for the effect of rotation. For the effec- 
tive radii we adopted the seeing-corrected values from Baggett, 
Baggett & Anderson (1998). In the few cases where this last 
compilation did not provide R e measurements (NGC 3982 and 
NGC 4800), we assumed a e = cr*. 

Figure 5 shows that, with one exception (NGC 4143), our de- 
rived upper limits on M B h lie within the scatter of the M B h - cr* 
relation, or above it. In this sense, the proposed M B h - cr* rela- 
tion passes this observational test for more than a dozen objects. 
Furthermore, except for the three galaxies with the lowest value 
of a*, our upper limits on Mbh exceed the values predicted by 
the M B h - cr* relation only by a modest amount (in average by 
a factor of ~ 4.6). 

In particular for the galaxies in our sample with actual M B h 
measurements from spatially resolved kinematics (S01), the 
present upper limits are consistent, within the errors, with the 
published M B h values. Moreover, when our +lcr upper lim- 
its are compared with other published M B h measurements (as 
compiled recently by Kormendy & Gebhardt 2001, Fig. 6), 
most appear to be quite close to the actual values for the black- 
hole mass, with the exception of the same galaxy lying below 
the M B h- cr* relation (NGC 4143) and the three objects in our 
sample with the lowest value of cr* (NGC 3982, NGC 4321, and 
NGC 4548). Figure 6 therefore shows that for the bulk of our 
sample with lOOkms" 1 < cr* < 200kms _1 our line-width mod- 
eling technique gives results statistically consistent with the val- 
ues obtained through other techniques. The current data do not 
provide any indication that the spiral and lenticular galaxies in 
our sample differ in their black-hole masses. Further, our upper 
limits on M B h do not seem to differ between barred or unbarred 
host galaxies, or as a function of the nuclear spectral classifica- 
tion. 

For any interpretation of upper limits, the basic sensitivity of 
the experiment is crucial. Our Mbh sensitivity limit does not 
only depend on the physical size of the resolution element and 
on the amount of stellar mass, but also on the spatial emissivity 
distribution of the gas. Indeed, we often find the spatial ex- 
tent of the ionized gas to be smaller than the dimension of the 
0"25 x 0."2 aperture. Considering the spatial emissivity of each 
of our sample galaxies, we can derive conservative sensitivity 
limits for a nearly face-on disk (cos; = 0.84) by first computing 
the predicted line widths arising in a purely stellar potential, 
and then by asking what values of Mbh need to be added in or- 
der to increase those line widths by their typical measurement 
error, conservatively ~10 km s" 1 (see §2.3). 
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FIG. 6. — Mbh upper limits versus Mbh determinations. In both panels the Mbh upper limits for our sample galaxies (symbols with downward arrows) are 
compared with the Mbh determinations that were deemed accurate in the recent compilation of Kormendy & Gebhardt (2001) (symbols without arrows). To be 
consistent with this compilation, for our sample objects with measured surface brightness fluctuation distance moduli (Tonry et al. 2001, for NGC 2787, NGC 3368, 
NGC 4143, NGC 4203, NGC 4459, and NGC 4548) the M B h upper limits have been rescaled assuming these distances. In the left panel the different symbols stand 
for the different techniques used to measure the Mbh, while in the right panel they represent the different types of host galaxies. The solid lines show the Mbh _ ct* 
relation as compiled by Gebhardt et al. (2000a). 



With a mean value of 3.9xl0 6 M Q , the derived sensitivity 
limits lie well below the inferred +la upper limits on Mbh (see 
Fig. 5 and Table 2), which can therefore be considered robust. 
We notice that by assigning a constant mass-to-light ratio to all 
Gaussian components of the luminosity density profile, includ- 
ing the ones with a < 0.5 pixel (see §3.3), the derived sensi- 
tivity limits increase only by 30%, to a mean value of 5.1 x 10 6 
M Q . 

Our +lcr upper limits for the whole sample correspond to 
Mbh values produced by nearly face-on disks (cos / = 0.84). As 
mentioned, we find one object, NGC 4143, for which the de- 
rived upper limit on M B h is below the M Bli - a* relation. For 
16 sample members, the expected number of disks with incli- 
nations more face-on than cos; = 0.84 is ~ 2-3. In order to 
reconcile the observed central line width of NGC 4143 with 
the Mbh _ c* relation within the Keplerian-disk framework, the 
disk would need to have an inclination angle of 27°. However, 
it should be noticed that NGC 4143 is among our sample one of 
the galaxies with a nearly unresolved spatial emissivity distri- 
bution (see Fig. 1), which could actually be more concentrated 
than the one adopted by us. Since this would lower the upper 
limit on M B h, we think that NGC 4143 represents an interesting 
candidate for future investigations. 

The situation is different for the three galaxies with the low- 
est values of a+ in our sample because their spatial flux profiles 
are resolved. Hence, we cannot explain their relatively high 
values of <j cen (<~ 100 kms" 1 ) within the context of a Keplerian 
disk in terms of gas orbiting in the vicinity of a ~ 2 x 10 6 M Q 
SMBH (as predicted by Mbh _ c* relation). Furthermore, as the 
derived sensitivity limits on M B h for these three galaxies are 
also around ~2x 10 6 M Q (see Table 2), we cannot expect the 
stellar mass contribution to help explaining their a cen values. 



Indeed, the line widths obtained from the stellar potential and a 
2 x 10 6 M Q SMBH are considerably smaller than the observed 
ones. The predicted line widths in this case are for these three 
galaxies in average <~ 40 and <~ 75 kms" 1 , in the nearly face-on 
(cos; = 0.84) and edge-on (cos; = 0.16) case, respectively. Al- 
ternatively, the observed central line widths might arise in all 
these three objects from highly inclined nuclear disks. This is 
not only unlikely but may also be insufficient, as in the case of 
NGC 3982, even when considering a perfectly edge-on nuclear 
disk. Therefore, in the case of the three less massive bulges, 
we may have indirect evidence that at least part of the observed 
line width is due to non-gravitational effects. 

5. DISCUSSION AND CONCLUSIONS 

We have demonstrated that with HSTs spatial resolution the 
integrated line widths of the central emission lines provide 
stringent and interesting constraints on the presence of SMBHs. 
The relative observational ease of this approach makes it poten- 
tially applicable to large galaxy samples, that would allow us 
to test the universal applicability of the emerging relations be- 
tween Mbh and galaxy properties. 

Our modeling, which was necessary to connect the observed 
o cen with the quantity of immediate interest, v c (R re f), was based 
on the assumption that the gas line width arises solely from or- 
bital motion within a randomly oriented disk around a putative 
SMBH. 

Reality is undoubtedly more complex, and we have consid- 
ered other potentially relevant effects, such as the stellar con- 
tribution to the total gravitational potential, and more simplisti- 
cally, hydrostatic support of the gas. The dynamical influence 
of outflows and magnetic fields could also be important. Ex- 
cept for fine-tuned circumstances, all these effects will provide 
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an additional contribution to the observed line width and the in- 
ferred upper limit on M Bli will be tighter. Hence, our adopted 
set of assumptions lead to conservative estimates. 

Comparison of our upper limits with direct Mbh determina- 
tions, either statistically (Fig. 6) or in a few cases individually 
(Fig. 5), showed that our l-a upper limits are generally near 
the actual value of M Bii . 

We have applied this analysis to a set of 16 galaxies whose 
sample selection was not biased toward particular M B h values. 
Remarkably, with one exception, our +la upper limits on Mbh 
closely parallel the M Bll -<j+ relation and suggest that for galax- 
ies with > 100 kms" 1 , SMBHs with exceptionally high M BH 
that violate the M B h - o* relation must be rare. By considerably 
broadening the range of host galaxies surveyed for SMBHs, our 
16 upper limits further support the emerging picture wherein 
the black-hole mass and the overall galaxy structure are closely 
linked. 

Even with a limited sample of 16 objects, we have been 
able to isolate a few cases worthy of further investigations. 



NGC 4143 stands out as the only object that falls below the 
Mbh - relation; we speculate that this may indicate that its 
nuclear disk is nearly face-on. Three low-ov galaxies (NGC 
3982, NGC 4321, and NGC 4548) seem to have M BH upper 
limits that lie systematically offset from other galaxies of low 
velocity dispersion in which the M B h was obtained by studying 
the stellar kinematics. This suggests that in low-mass bulges 
non-gravitational forces can considerably affect the gas motions 
in the central 10 pc. 
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Table 1 

Basic Parameters of the Sample Galaxies 



Galaxy 


Hubble Type 


Bj 


Spectral Class 


D 




Ref. 


<J e 


&cen 


Obs. Date 






(mag) 




Mpc 


(km s" 1 ) 




(km s" 1 ) 


(kms" 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 2787 


SB0 + 


11.82 


L1.9 


13.0 


210±23 


1 


185 ±20 


215.1 ±4.3 


05 Dec. 1998 


NGC 3351 


SBb 


10.53 


H 


8.1 


101 ±16 


2 


93 ±15 


47.2 ±1.5 


25 Dec. 1998 


NGC 3368 


SABab 


10.11 


L2 


8.1 


135±10 


3 


114±8 


101.5 ±3.1 


31 Oct. 1998 


NGC 3982 


SABb: 




S1.9 


17.0 


78 ±2 


4 


78±2 


136.7±3.5 


11 Apr. 1998 


NGC 3992 


SBbc 


10.60 


T2: 


17.0 


140 ±20 


4 


119± 17 


109.5 ±3.1 


19 Feb. 1999 


NGC 4143 


SAB0° 


11.65 


L1.9 


17.0 


270 ±12 


5 


271 ±12 


226.3 ±2.1 


20 Mar. 1999 


NGC 4203 


SABO": 


11.80 


L1.9 


9.7 


124±16 


1 


110± 14 


148.9±5.2 


18 Apr. 1999 


NGC 4321 


SABbc 


10.05 


T2 


16.8 


83 ±12 


6 


74± 11 


85.7±1.7 


23 Apr. 1999 


NGC 4450 


Sab 


10.90 


L1.9 


16.8 


130± 17 


2 


121 ±16 


162.4 ±1.7 


31 Jan. 1999 


NGC 4459 


S0 + 


11.32 


T2: 


16.8 


189±21 


1 


167± 18 


193.1 ±5.2 


23 Apr. 1999 


NGC 4477 


SBO:? 


11.38 


S2 


16.8 


156± 12 


7 


134±10 


128.9±2.2 


23 Apr. 1999 


NGC 4501 


Sb 


10.36 


S2 


16.8 


151 ± 17 


8 


136± 15 


110.8±1.8 


26 Apr. 1999 


NGC 4548 


SBb 


10.96 


L2 


16.8 


82±9 


9 


71 ±8 


81.2± 1.8 


26 Apr. 1999 


NGC 4596 


SB0 + 


11.35 


L2:: 


16.8 


154 ±5 


10 


136±4 


142.0 ±8.7 


20 Dec. 1998 


NGC 4698 


Sab 


11.46 


S2 


16.8 


134±6 


9 


116±5 


101.9±2.2 


24 Nov. 1997 


NGC 4800 


Sb 


12.30 


H 


15.2 


112±2 


4 


112±2 


71.9±7.0 


03 Mar. 1999 



Note. — Col. (1): Galaxy name. Col. (2): Hubble type from de Vaucouleurs et al. (1991). Col. (3): Total apparent B 
magnitude from de Vaucouleurs et al. (1991). Col. (4): Nuclear spectral class from Ho et al. (1997), where H = H II nucleus, 
L = LINER, S = Seyfert, T = "transition object" (LINER/H II), 1 = type 1,2 = type 2, and a fractional number between 1 and 
2 denotes various intermediate types; uncertain and highly uncertain classifications are followed by a single and double colon, 
respectively. Col. (5): Distance from Tully (1988), who assumes Ho = 75 kms" 1 Mpc" 1 . Col. (6): Ground-based central stellar 
velocity dispersion cr*. Col. (7): Reference for cr*. Col. (8) Stellar velocity dispersion within one effective radius a e , derived from 
the value following the extrapolation of J0rgensen, Franx, & Kjaergaard (1995). Col. (9): Measured gas velocity dispersion 
within the central 0"25 x 0."2. Col. (10): UT observation date. 

References. — (1) Dalle Ore et al. 1991; (2) Whitmore, Schechter, & Kirshner 1979; (3) Heraudeau et al. 1999; (4) this paper 
(§2.2); (5) Di Nella et al. 1995; (6) Schechter 1983; (7) Jarvis et al. 1988; (8) Heraudeau & Simien 1998; (9) Corsini et al. 1999; 
(10) Kent 1990. 
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Table 2 

Results of the Rotating Disk Model 



Galaxy Name 


Pflux 


&flux,corr 


Mbh($ = $bh) 


M bh ($ = $bh + $*) 


M B H(Sensitivity Limit) 




(arcsec) 


(arcsec) 


(M ) 


(M ) 


(M ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 2787 


0.056 


0.053 


1.9xl0 8 


1.8xl0 8 


3.8xl0 6 


NGC 3351 


0.132 


0.153 


9.7xl0 6 


8.0xl0 6 


2.1 xlO 6 


NGC 3368 


0.086 


0.095 


3.8xl0 7 


2.7 xlO 7 


4.5 xlO 6 


NGC 3982 


0.051 


0.046 


8.0xl0 7 


7.5 xlO 7 


3.3xl0 6 


NGC 3992 


0.051 


0.046 


5.7xl0 7 


5.3xl0 7 


3.0xl0 6 


NGC 4143 


0.035 


0.024 


1.4xl0 8 


1.4xl0 8 


1.4xl0 6 


NGC 4203 


0.030 


0.015 


2.3xl0 7 


2.3 xlO 7 


4.3 xlO 5 


NGC 4321 


0.046 


0.040 


2.7xl0 7 


2.5 xlO 7 


2.3 xlO 6 


NGC 4450 


0.046 


0.040 


1.1x10 s 


l.lxlO 8 


3.1 xlO 6 


NGC 4459 


0.041 


0.033 


1.3xl0 8 


1.3xl0 8 


2.9xl0 6 


NGC 4477 


0.056 


0.053 


8.7xl0 7 


7.8 xlO 7 


5.1 xlO 6 


NGC 4501 


0.081 


0.088 


9.0xl0 7 


7.4xl0 7 


7.6xl0 6 


NGC 4548 


0.066 


0.067 


3.3xl0 7 


2.8 xlO 7 


3.8xl0 6 


NGC 4596 


0.056 


0.053 


l.lxlO 8 


9.4xl0 7 


5.6xl0 6 


NGC 4698 


0.091 


0.100 


8.0xl0 7 


7.1 xlO 7 


6.0xl0 6 


NGC 4800 


0.076 


0.081 


3.2xl0 7 


2.0xl0 7 


6.1 xlO 6 



Note. — Col. (1): Galaxy name. Col. (2): Derived first guess of the intrinsic flux extent(see §3.2). Col. 
(3): Adopted intrinsic flux extent(see §3.2). Col (4): Upper limit onM B H (+lc) obtained with a rotating disk 
model ignoring the stellar gravitational potential (§3.2). Col (5): Upper limit on M B h (+1ct) obtained with a 
rotating disk model including the stellar gravitational potential (§3.2). Col. (6) Sensitivity limit to M Bll - 



